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A TIMER INVENTORY BASED UPON MANUAL AND AUTOMATED ANALYSIS OF 
ERTS-1 AND SUPPORTING AIRCRAFT DATA USING MULTISTAGE PROBABILITY 
SAMPLING 

James D. Nichols, Mike Gialdini and Sipi Jaakkola, Center for Remote Sensing Research, 260 Space Sciences 
Laboratory, University of California, Berkeley, California 


ABSTRACT 

A quasi-operational study demonstrating that a timber inventory 
based on manual and automated analysis of ERTS-1, supporting air- 
craft data and ground data was made using multistage sampling 
techniques . 

The inventory proved to be a timely, cost-effective alternative 
to conventional timber inventory techniques. The timber volume 
on the Quincy Ranger District of the Plumas National Forest was 
estimated to be 2.44 billion board feet with a sampling error of 
8.2 percent. Costs per acre for the inventory procedure at 1.1 
cent/acre compared favorably with the costs of a conventional 
inventory at 25 cents/acre. A point-by-point comparison of 
CALSCAN-classified ERTS data with human-interpreted low altitude 
photo plots indicated no significant differences in the overall 
classification accuracies. 


SECTION I: INTRODUCTION AND PROCEDURES 

In order to test the usefulness of ERTS imagery for wildland re- 
source inventories , a timber inventory was performed in which the ERTS 
imagery acted as the first stage of a multistage sampling design. The 
objective of the inventory was to estimate the standing volume of mer- 
chantable timber within the Quincy Ranger District (215,000 acres) of 
the Plumas National Forest in California . Secondary ob j ectives of the 
inventory were: (1) to test the operational efficiency of the sampling 

procedures of the multistage sampling design; (2) to test the effective- 
ness of the CALSCAN classifier on the ERTS data; (3) to determine the 
value of ERTS data and the aircraft data in reducing the sample error; 
and (4) to compare the costs of this timber inventory with other inven- 
tories that utilize conventional procedures. 

A three stage sampling design was tested in which "timber volume" 
was the variable estimated. At each stage timber volume estimates were 
made from sampling units whose probability of selection in the sample 
was proportional to the predicted volume. Timber volume estimates were 
made from three stages: (1) the first stage involved automatic classi- 

fication of the timberland on the ERTS data tapes into four timber voi- 
um classes. Within the classified area subsamples were selected 
(called primary sampling units; PSU) from which a more refined estimate 
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of timber volume could be made in the second stage; (2) the second 
stage involved the acquisition of low altitude photograph of selected 
primary sampling units to choose photo plots based on a second timber 
estimate made by comparison with photo-volume tables; (3) the third 
stage involved selecting individual trees within selected sample photo 
plots by photo measurement of all merchantable trees. The selected 
trees were then precisely measured for volume on the ground and these 
volume measurements in turn were expanded through the various stages 
of the sample design to estimate total timber volume over the national 
forest land within the Quincy Ranger District. 

The statistical procedures for expanding the timber volume esti- 
mates through the various stages of the timber inventory are discussed 
in Section II of this study. The procedures for selecting sampling 
sites and estimating timber volumes from ERTS and aircraft imagery at 
each of the stages of the sample design are described in subsequent 
sections in order to demonstrate how the ERTS and supporting aircraft 
data were necessary components in performing the timber inventory. 

Stage I. CALSCAN Classification of ERTS Data and Primary Sample Unit 
Selection 
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ERTS-l data tapes of the Quincy Ranger District, Plumas National 
Forest, were classified on the CRSR interactive human-computer system 
using a CALSCAN point-by-point classification routine. The coordinates 
of the Ranger District boundary and those of non-national forest land 
within the District were identified on the tapes so that only those 
picture elements associated with national forest land were classified 
and incorporated into the inventory. This procedure considerably re- 
duces the costs of classification. 

The classification was based on four timber volume classes, namely 
(1) non forest; (2) forest sites containing less than 10,000 Bd. ft./ac.; 
(3) forest sites containing 10,000-20,000 Bd. ft./ac.; and (4) forest 
sites containing more than 20,000 Bd. ft./ac. The classifier was trained 
to recognize each of the four timber volume classes based upon photo- 
interpreter selection of 33 training cells whose range of timber volumes 
was based on crown closure and average crown diameter. The training 
cells were selected from interpretation of high altitude color-infrared 
photography (scale 1/120,000). Each of the training cells was located 
oi the ERTS imagery and digitized. Point-by-point classification of 
all ERTS data points within the Quincy Ranger District proceeded by 
matching each data point (picture element) with the training cells. The 
results were grouped into the four timber volume classes. The accuracy 
achieved by automatic classification of ERTS tapes (Tables 1 and 2) 
demonstrates the efficiency which can be attained in the timber in- 
ventory through analysis of ERTS data tapes in the initial stages of the 
sample design. 
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The classified ERTS data (CALSCAN classification) of the Quincy 
Ranger District was divided into rectangular sampling units (called 
primary sampling units). Each unit measured 1325 ft. wide by one and 
a half miles long. The size of these sampling units was based upon a 
practical area which could be photographed in a single flight line by 
a light aircraft using a 35 mm camera system, the ability of the ground 
crew to complete the ground work for a flight line in one day, and the 
variation between sampling units. For each primary sampling unit, the 
following information was computed. 

(1) The number of points in each volume class (within the unit). 

(2) The weighted total volume for each volume class.* 

(3) The sum of the weighted totals for all classes. 

(4) A cumulative sum of the weighted totals. 

(5) The mean volume for all sampling units. 

(6) The variance of the sampling units. 

Based upon the information either estimated or computed for each 
primary sampling unit (PSU) , four units were selected for further sam- 
pling in the timber inventory. The four units were selected with prob- 
ability of selection proportional to their estimated volumes. The 
location of the four selected PSU's was transferred from the ERTS clas- 
sified images to the color-infrared high altitude aerial photography 
(scale 1/120,000) to facilitate locating them accurately from the 
air when they are photographed from a lower altitude as part of the 
second stage of the timber inventory. 

Stage II. Volume Estimation on Low Altitude Photography 

Two 35 mm cameras were used to obtain low altitude photography of 
the selected primary sampling units at two different scales. A 24 mm 
focal length wide-angle lens was used to acquire complete coverage of 
each sampling unit at an approximate scale of 1/7500, and a 200 mm focal 
length lens obtained large scale stereo triplicates, scale approximately 
1/1000, from which to make precise photo estimates of timber volume. 

The camera with the telephoto lens was equipped with a motorized film 
drive which enabled each stereo triplicate to be taken within one second 
at five second intervals , while the camera with the wide-angle lens was 
operated manually to obtain single frames at five second intervals. The 
photo coverage for each PSU consisted of 10 stereo triplicates and 10 
wide-angle photographs. 

The wide-angle photos of each primary sampling unit were formed into 
a mosaic to show its full area. The center of the middle photo for each 
stereo triplicate was used as a plot center and was located and marked 
on the mosaic. Each plot center was also located on a topographic map 
and its elevation was determined. 


*The weighted total volume was determined by multiplying the number 
of points in each volume by the assigned weight for that class : 
non-timber = 0; less than 10,000 Bd. ft./ac. = 1; 10,000-20,000 
Bd. ft./ac. = 2; 20,000 Bd. ft./ac. and above = 3. 
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Discriminant Analysis 



NF = Non forest 

T1 = Less than 10,000 
board feet per acre 

T2 = 10,000 to 20,000 
board feet per acre 

T3 = Greater than 20,000 
board feet per acre 


Table 1 

The accuracy obtained in a comparison of classified ERTS data with 
photo interpreter classification from a photo scale of 1/1000. The dia- 
gonal of the matrix shows the percentage correct classification when 
comparing the discriminate analysis of ERTS data with the large scale 
photographs. The values not along the diagonal represent the percentage 
of points which were classified into adjacent timber volume classes. 

The values at the bottom of each column indicate the percentage of 
points classified by discriminate analysis of ERTS data which were great- 
er than (+) or less than (-) the points classified by interpreters on 
large scale aerial photographs. 
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Discriminant Analysis of ERTS-1 data 

Com Soil/R BR H.W Wot Conv Gross CONI s Commercial conifer consisting of 
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The scale of each photo plot was determined, and a .4 acre circular 
plot was drawn about the photo plot center. Timber volume in each .4 
acre photo plot was estimated by referring to photo-volume tables based 
upon interpretation of % crown closure and measurement of average stand 
height using a parallax bar (Chapman, 1965) . Within each primary sam- 
pling unit, two out of the ten possible photo plots were chosen with 
probability of selection proportional to the estimated volume. 

Stage III. Selection of Trees for Precise Ground Measurement of Timber 
Volume 

In the third stage, all trees of merchantable size within each 
selected photo plot were pin-pricked and numbered. For each of these 
trees, the average crown diameter was determined based on the longest 
and shortest dimension of their crowns. After adjustments for scale, 
the average crown diameter value was cubed (raised to the third power) 
to be used as a relative measure of individual tree volume for the 
third stage volume estimation. Four trees were selected from the popu- 
lation of merchantable trees found within each photo plot to be measured 
by a dendrometer on the ground. (Selection again was based upon proba- 
bility proportional to the estimated volume of each tree.) 

A two-man crew went into the Quincy Ranger District with a Barr 
and Stroud optical dendrometer to measure the selected trees. The 
large scale (low altitude) photographs were used to locate the photo 
plot centers as well as the trees within the plots to be measured. In 
addition to the dendrometer measurements, an easily recognizable feature 
on the ground near the plot center was measured in order to get a more 
accurate estimation of the photo scale of the plot. The dendrometer 
measurements were brought back from the Quincy Ranger District and en- 
tered into a computer program that calculated volumes for the individual 
trees. The tree volumes were then expanded through each stage of the 
sample design to estimate total volume on the District, consistent with 
the statistical methods for variable probability sampling (see Section 
II). 


SECTION II. STATISTICAL METHODS FOR TIMBER VOLUME ESTIMATION 

Timber volume predictions were made from three stages of the timber 
inventory of the Quincy Ranger District for the purpose of selecting 
sample plots whose probability of selection was proportional to the 
volume predictions. Thus, variable probability sampling methods were 
used to estimate the total volume in this timber inventory. Three va- 
riables proportional to timber volume were used in generating the se- 
lection probabilities: (1) "Volume" estimate of the ERTS picture ele- 

ment based on the spectral signatures on four bands and subsequent 
training and classification; (2) volume estimate of a plot on a 1:1000 
scale color print based on the photo-volume tables (Chapman, 1965) ; 
and (3) volume estimate on a large scale photo based on a rough volume 
estimate - crown diameter cubed. 
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Using a scheme where probability of selection is proportional to 
the estimated volume, the effort is focused on the areas of higher tim- 
ber volume and adds to the overall cost-efficiency. The ability to 
list the populations at each stage prompted the selection of list sam- 
pling as the variable probability sampling scheme. 

Method of Estimation 

The method of estimation was based on "unequal expansion" as im- 
plied by the probability scheme discussed above. At each of the three 
stages, the probability proportional to estimated size (p^) was obtained 
by listing the volume estimates of the sampling units (x^) , and dividing 
them by the total of volume estimates/ n \ : 

(,?, ' 


n 

L x i 

i=l 


A sample of a chosen size was then drawn by applying random inte- 
gers from 1 to n x^ and observing the probability interval and the 
corresponding ^ sampling unit which contains the randomly selected 
integers. i=l 

In the first and second stages the timber volumes of the selected 
sampling units were estimated by subsequent sampling, whereas in the 
third stage the volume was carefully measured by a precision dendrometer. 
The entire three stage estimation procedure was as follows: 

Stage I: A sample of n^ out of the PSU's was drawn from stratum 

h with probability proportional to estimated size (ppes) . The estimate 
of the total volume then becomes: 



where: L = total number of strata 

P = selection of probability of the 

1 ith PSU in the hth stratum 


151 



= total volume of the 1th PSU in 
the hth stratum (remains to be 
estimated by subsequent stages) . 

Stage II : To estimate the total volume (y, .) of the ith PSU, a sam- 

ple of n,. out of the secondary sampling units (.4 acre plots) is 
drawn with probability proportional to estimated volume. This would 
give: 


y 


hi 


Hhi 

S Jk 

J-X Phi3 


However, in order to include area expansion from circular sample 
plots to the full PSU, plus stratify the second stage plots into four 
volume strata, the estimator becomes: 




\ir 

^ir 


_1 

n hir 


Hfiir 

E 


i=i 


A 

Y hirj 

P hirj 


where: r = 1,2,...R refers to the CALSCAN 

volume strata 


P^ ir = selection probability of the rth 
volume stratum of the ith PSU in 
the h th stratum 

A = area (indices as above) 

a = sample area (indices as above) 

n = sample size (indices as above) 


P, . , = selection probability of the jth 
J plot of the r th volume stratum, of 
the ith PSU in the hth stratum 



plot volume (to be estimated by 
Stage III) 


Stage III : To estimate the total volume of the jth plot , a sample 
of nj^ r . out of the N, . . tertiary sampling units (trees) is drawn with 
ppes. Then ... J 
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A 

y hlrj 


1 

“hirj 


n hirj 

z 

k=l 


Y hirjk 

P hirjk 


where: P h±rjk 


the selection probability of the 
kth sample tree of the jth plot of 
the rth volume stratum of the ith 
PSU of the hth stratum. 


Y hirjk 


the dendrometer measured volume of 
the kth sample tree of the jtfi plot 
of the r th volume stratum of the 
ith PSU of the hth stratum. 


Combining the various stages above, the entire estimator becomes: 


L 

Z 


i_ 

% 



1 \ir 1_ 

P hir ^ir ^ir 



1 1 
P hirj ^irj 



k=l 


y hirjk 

P hirjk 


Variance of the Estimator 


In multistage sampling, when the number of first stage units is 
large, most of the variability in the population is due to the first 
stage. Therefore, it suffices to consider only the first stage values 
(here y, ,) to estimate the population variance and, consequently, the 
variance of the estimator (Durbin, 1953, p. 262; Kendall and Stuart, 
1967, vol. 3, p. 200; Langley, 1971, p. 131). 


Thus for the first stage our stratified sampling estimator becomes 
(Cochran, 1963, p. 260): 


A 

V 


L 


z 

h=l 


1 ^ 

nh 



y 


hi 


P 


hi 


Its variance is: 


Var 


A 


(V) 



which has an unbiased estimator 


Var (V) 


L 


£ 

h=l 


n h (n h " 1} 




For proportional allocation, n^ = n (N^/N) and 


Var (^) 
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£ 

h=i 


_N 

n N, 
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Var (V) 




h=l n N h (nN h - 




The last equation is an unbiased estimator of Var (V) and can be 
used for the estimation of the sampling error of the inventory. 


SECTION III. RESULTS AND CONCLUSIONS 

The total volume of timber in the Quincy Ranger District of the 
Plumas National Forest was estimated to be 407 million cubic feet (ap- 
proximately 2.44 billion board feet) based on eight selected photo 
plots located within four primary sampling units. The sampling error 
associated with this estimate was 8.2% which falls below the expected 
sampling error of 20% for the Quincy Ranger District. This indicates 
that the true volume of merchantable trees in the Quincy Ranger District 
will fall into the interval 2.11 - 2.77 billion board feet with 80% 
probability. 

There were only 31 trees total measured by an optical dendrometer 
on the ground at the eight plots (32 trees should have been measured 
but one plot out of the eight contained only three merchantable trees 
which could be measured). The field work required one week's time from 
a two-man crew, and the total area of the ground plots measured was 3.2 
acres, representing a sampling fraction of approximately 1/67,000. 

Table 1 lists the costs of the timber inventory on 215,000 acres 
using the multistage sampling design as an operational system. The 
expected costs of an inventory using the above procedures, on the en- 
tire Plumas National Forest (1,161,554 acres) would be approximately 
$15,000 and would take 5 months to complete. In comparison, the U. S. 
Forest Service in 1970 completed an inventory of the Plumas National 
Forest using the 10 point system, a conventional procedure, at a cost 
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of $300,000 and took two years to complete. Some deficiencies of that 
inventory were: (1) that it did not provide in-place mapping of timber 

volumes and (2) it used interpolation to arrive at the volumes for each 
ranger district which tends to reduce the accuracy of the volume esti- 
mates for each district, as compared to the estimate for the entire fo- 
rest. The CRSR system would inventory each ranger district as a sepa- 
rate unit, thereby achieving more precise estimates and could also pro- 
vide in-place mapping of volume. 

The Forest Service inventoried the entire El Dorado National Forest 
(858,496 acres) using WRIS (Wildland Resource Inventory System) at a 
cost of approximately $31,300. While WRIS closely resembles the CRSR 
system in design, the results again lacked in-place mapping of volume 
as well as not treating each ranger district as a stratum. In terms of 
inventory costs per acre for each system, the CRSR system appears to be 
best (10-point system = 25c/acre, WRIS = 3.6c/acre, and CRSR = 1.1 0/ 
acre) while still providing more information with little or no loss in 
precision. It should also be noted that the multistage design used 
here becomes more efficient as the total area considered increases. 

For example, if the entire State of California were inventoried, the 
per acre costs would be . 175<:. 
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TABLE 3 

COST ESTIMATE-TIMBER INVENTORY 
BASED ON 215,000 ACRES - QUINCY RANGER DISTRICT 


Computer Processing 


CALSCAN Training and Photo Interpretation 

$ 120 

CALSCAN Classification 

210 

Statistical Break-up 

40 

CALSCAN Statistics 

12 

Aerial Photography 


Aircraft and Crew 

$ 210 

Film and Processing 

90 

Supplies and Expenses 


Travel 

$ 350 

Miscellaneous 

50 

Personnel 


1 Project Scientist 

2 mos. F.T.E. @ $1100/mo. 

$2,200 

1 Software Consultant 

.5 mo. F.T.E. @ $1300/mo. 

550 

1 Statistician 

.5 mo. F.T.E. @ $1100/mo. 

550 

2 Lab Assistants 

1 mo. F.T.E. Each @ $620/mo. 

Salary Subtotal 

1,240 

$4,640 

Associated Overhead 

TOTAL 


$ 382 


$ 300 


$ 400 


$968 


$5,568 


$6,640 
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SUMMARY 


The preliminary results of the timber inventory of the Quincy 
Ranger District indicate that the procedures employed in the multi- 
stage sampling design are valid and substantially reduce both the 
costs and the amount of time required to perform a timber inventory 
for a large area. This study demonstrates the value of ERTS data 
for accurately correlating picture elements with timber volume esti- 
mates as a fundamental first step in selecting primary sampling units 
in the first stage of the inventory. The inventory procedures uti- 
lized here will be applied to two additional districts of the Plumas 
National Forest in an effort to estimate the total timber volume on 
the forest to a desired sampling error of 10%. 

Future activities will include the investigation of the ability 
of the CALSCAN classifier to correlate timber condition classes with 
picture elements to provide information for use in forest simulators. 
This would enable forest managers to project long-range results of 
management practices to determine their effectiveness. Also, similar 
inventory procedures will be developed for application to other vege- 
tation types to provide timely, cost-effective, and useful information 
for resource management planning, such as fuel mapping for fire control 
and management and rangeland-condition mapping for carrying-capacity 
determinations . 
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